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Under natural conditions, excitation of biological molecules, which display non-unitary open 
system dynamics, occurs via incoherent processes such as temperature changes or irradiation by 
sunlight /moonlight. The dynamics of such processes is explored analytically in a non-Markovian 
generic model. Specifically, a system S in equilibrium with a thermal bath TB is subjected to an 
external incoherent perturbation BB (such as sunlight) or another thermal bath TB', which induces 
time evolution in (S+TB). Particular focus is on (i) the extent to which the resultant dynamics is 
coherent, and (ii) the role of "stationary coherences", established in the (S+TB) equilibration, in 
the response to the second incoherent perturbation. Results for systems with parameters analogous 
to those in light harvesting molecules in photosynthesis show that the resultant dynamical behavior 
is incoherent beyond a very short response to the turn-on of the perturbation. 
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I. INTRODUCTION 

The dynamics of open quantum systems (e.g., molecules 
of interest in contact with an environment) in the presence 
of an external perturbation is of great current interest 
and applicability. Two general circumstances can be envi- 
sioned, one where the external perturbation is "designed" , 
the second where it occurs naturally. Within the frame- 
work of electromagnetic perturbations, popular examples 
of the former include pulsed laser studies of molecular 
dynamics in the laboratory [1, 2], laser based scenarios 
for quantum control and quantum information, various 
spectroscopies, etc. These artificial light sources are often 
characterized by pulses of short temporal duration and 
significant coherence times. 

By contrast, natural light induced processes such as vi- 
sion or photosynthesis are induced by virtually stationary 
chaotic blackbody light sources such as the sun, which 
have significantly different properties [3-7]. Similarly, nat- 
ural processes, like ion transport through membranes[8] 
may be induced by temperature changes. Both perturba- 
tions, temperature change and excitation by natural light, 
are fully incoherent. That is, they are both associated, as 
outlined later below, with a perturbation described by a 
density matrix that is a mixture of stationary states. 

A proper treatment of the systems described above 
requires that the system S of interest be first equilibrated 
with a background thermal bath TB, after which it is 
subjected to a second perturbation comprising a second 
thermal bath (denoted TB') or blackbody radiation (de- 
noted BB). Our focus here, in this two step process, is on 
the role of coherences in the dynamics. 

One motivation for this work lies in recent studies of 
coherent quantum dynamics in model photosynthetic light 
harvesting systems [1, 9] and in vision[2, 4, 5]. Coherent 
dynamics of this kind is observed in experiments in which 
the system is excited with coherent laser light, and the 
timescales for the decay of the coherence, generated by the 



system interacting with its environment, is then measured. 
Principal amongst the observations on molecules involved 
in photosynthesis is that the coherent dynamics associ- 
ated with electronic energy transfer persists, at room 
temperature, on considerably longer time scales (e.g., 400 
fs - 2000 fs) than is expected from earlier results on other 
systems[10-13], and from theories of decoherence[14]. A 
number of computational results have been obtained, and 
theories advanced, as to why such longevity occurs. For 
example, we have identified a number of physical condi- 
tions under which such long lived electronic coherence 
persists[15, 16]. 

The relevance of these results for realistic photosyn- 
thetic systems, or for the operation of devices that mimic 
photosynthesis, however, depends heavily on the relation- 
ship between the dynamics observed in the experiments, 
which use pulsed coherent light, and dynamics under nat- 
ural light, such as that from the sun. Both recent [7] as 
well as earlier studies[3] on this relationship for isolated 
molecules (i.e. molecules that are not in contact with an 
environment) show that whereas pulsed coherent light in- 
duces dynamics in the molecule, natural stationary chaotic 
light does not. Rather, irradiating an isolated molecule 
with such natural light over long natural time periods 
creates a stationary mixture of molecular eigenstates. 

Natural light absorbing molecules of interest in, e.g., 
photosynthesis and vision are not, however, isolated. 
Rather they are in contact with an external environ- 
ment through which decoherence and relaxation occur. 
Hence, there is a need to understand the response of 
a system+bath that is subject to an incoherent excita- 
tion. A study of Retinal excitation[5] under these con- 
ditions confirmed that stationary eigenstates resulted in 
this open system as well. An alternate study on photo- 
synthetic systems [6] simultaneously in contact with both 
(TB+BB) noted the appearance of off-diagonal elements 
of the subsystem density matrix S and regarded them 
as "stationary coherences" , a concept elucidated further 
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below. The relation of these stationary coherences to 
dynamics was, however, unclear. Hence, clarifying the 
situation regarding open systems under incoherent per- 
turbation is, therefore, well motivated. 

As noted above and as emphasized, in this work it is 
important to recognize that the natural processes noted 
above occur in two distinct steps. That is, in the first step 
the subsystem S thermally equilibrates with its surround- 
ing bath, yielding a result that we denote as (S+TB). 
In the second step this equilibrated system is placed in 
contact with a source of incoherent light (BB) or another 
thermal bath TB' with which the system interacts, giving 
an (S+TB)+BB [or (S+TB) +TB'] system. As described 
in detail below in Sec. II, the first generates "stationary 
coherences" . The nature of the subsequent relaxation dy- 
namics, the role of the stationary coherences arising from 
the equilibration of S with TB in the subsequent dynamics 
of (S+TB)+BB [or of (S+TB)+TB'] are described below. 
Of particular interest are the nature and time scales over 
which the dynamical and equilibration process occurs, 
and the role played by the stationary coherences. 

The existence of substantial stationary coherences re- 
quires strong system-environment coupling (whether dis- 
sipative or not) [17, 18], and as we will show, are largest 
in the low temperature regime. This regime of strong 
coupling and low temperature is relevant in the context 
of, e.g., electronic energy transfer[15, 16]. Being in this 
regime prevents the application of simplifying approxima- 
tions to the evolution of the density matrix, such as the 
Markovian approximation[19] or secular approximation. 
That is why a thorough analysis of these off-diagonal 
terms and their influence on the dynamics have been elu- 
sive. Here we show, using an analytic model that these 
off-diagonal terms can contribute to the dynamics of the 
populations only when they enter as initial correlations 
between the system and the thermal bath, and that in 
realistic systems these contributions are small. 

Note for clarity below, the specific characteristics of 
the natural process that are of interest. That is, we are 
focusing on molecular systems that are, as in natural pho- 
tosynthetic processes, irradiated by blackbody radiation 
(or another thermal source) for time scales that are far 
far longer (e.g. hours) than the inverse of the molecular 
energy level spacing, which would define the time scale 
for coherent molecular dynamics. Results of this study 
are found to extend the result previously obtained for the 
isolated system to the open system. That is, "natural 
incoherent chaotic light" , e.g. from the sun or the moon, 
is shown to be incapable of inducing coherent dynamics 
in either isolated or open systems. 

Note that the dynamics examined here takes place on 
a single electronic surface. The extension of this analytic 
influence functional approach to non-adiabatic processes 
is in progress. However, the general conclusions obtained 
here are expected to hold in the case of excitation from 
one electronic surface to another. 

The paper is organized as follows: Sec. II introduces 
features of the first thermalization step, i.e. the process of 



S+TB relaxing, with a focus on the concept of "stationary 
coherences" . Section III provides the Hamiltonian under 
consideration and associated computational results for the 
(S+TB) relaxation. The case of (S+TB) interacting with 
a second thermal bath is contained in Sec. IV, whereas 
the case where the perturbation is incoherent light is 
discussed in Sec. V. A short discussion and summary is 
provided in Sec. VI. 



II. EQUILIBRATION OF (S+TB) AND 
"STATIONARY COHERENCES" 

The essence of open quantum systems lies in the in- 
teraction between a part of interest, "the system", with 
a component which is not of interest "the environment" 
or "bath". Consider then the two components, a system 
Hg with eigenbasis {rii} defined on the Hilbert space 
Hs and the thermal bath (in this section), described by 
the Hamiltonian Htb with eigenbasis {-/Vj} defined on 
the Hilbert space "Htb • In the absence of an interaction 
between them, the total Hilbert space T-L = Hs <8> "Htb 
will be diagonal in the basis {rii ® Ni}. However, if the 
systems interact via a coupling term ifsT, then % is no 
longer diagonal in this basis. To obtain information on the 
system we would then trace over the degrees of freedom 
of the bath. It is then clear that the remaining sector 
"Hs is no longer diagonal in the {rii} basis. Rather, the 
new effective eigenbasis should "know about the ther- 
mal bath" ; therefore the resultant effective system basis 
should be a function of the coupling to the bath and of 
the temperature. 

In general, one is not interested in studying the dynam- 
ics in terms of this coupling-and-temperature dependent 
basis, but rather in describing the dynamics in the sys- 
tem eigenbasis of H^. Viewed in this basis, the coupling 
to a bath generates off-diagonal elements in the system 
density matrix. Once the (system+bath) have relaxed 
to equilibrium these off-diagonal elements do not change 
with time, and are not associated with dynamics. Rather, 
they can be termed "stationary coherences". They are 
just a manifestation of our focus on the system, rather 
than on the coupled system+bath. 

An alternate perspective on these stationary coherences 
is that strong coupling to the environment causes over- 
lap of homogeneous line shapes associated with different 
energy states of system S. 

Given that the equilibrated S+TB composite is the 
natural state in which one finds the system S in the cases 
of interest, the relevant questions here are then: (i) what 
is the dynamics when this equilibrated (system+bath) is 
exposed to an incoherent perturbation, and (i) what is 
the role of the stationary coherences in this subsequent 
dynamics? For example, it has been suggested[6], but not 
explored quantitatively, that such terms are capable of 
inducing coherent dynamics when the equilibrated S+TB 
is subject to an additional perturbation. 

Below we show how to properly address the incoherent 
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excitation of a system S in a thermal environment by, in 
contrast with previous work [4-6], considering the exci- 
tation in two steps. In particular, using an analytically 
soluble model valid in the whole range of parameters, 
we study the role of environment-generated coherences 
during: (i) thermalization of an initially isolated cen- 
tral system S in contact with a thermal bath TB, and 
(m) the subsequent dynamics induced by the presence 
of a second thermal bath TB' at different temperature 
or by blackbody radiation BB. In the absence of initial 
coherences in the system eigenbasis, we show in situation 
(z) that, although time dependent off-diagonal elements 
could be detected, they are unrelated to the dynamics of 
the diagonal terms. In situation (u), coherences having 
been initially generated in S by TB contribute naturally 
to the dynamics of the populations of S toward the new 
equilibrium state. 



III. COHERENCES AND DYNAMICS TOWARD 
THERMALIZATION 

The process we described above is not a particular 
feature of light-harvesting systems, but is a generic fea- 
ture of open quantum systems and is therefore ubiqui- 
tous. In order to appreciate this in detail, and to ex- 
tract the relevant features of the underlying physical 
situation, we consider a model: a harmonic oscillator 
immersed in a dissipative environment TB [20-22]. Al- 
though an idealization, it encompasses a reasonable de- 
scription of a wide variety of objects in nature such as 
low energy vibrational molecular modes, in addition to 
artificial ones like nanomechanical oscillators, optical and 
microwave cavities, and movable mirrors. [23]. The Hamil- 
tonian of the system+environment can be written as 
H = Hg + Htb + -ffsTj where Hg is the Hamiltonian of 
the unperturbed oscillator, Htb is the Hamiltonian of 
the thermal bath TB, and H$t describes the interaction 
of the system with TB. In particular, we choose 
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and Hs — + ^j^Qx, where m and uiq are the natu- 

ral mass and frequency of the oscillator while p x and q x 
denote its (the thermal bath mode j) canonical conjugate 
momentum and position coordinates. Momenta and co- 
ordinates of bath mode j are denoted pj and qj and Cj 
denotes the coupling term to the j-th mode. In Fig. 1, 
we have depicted the physical situation. 

The evolution of the system density-matrix can be 
analytically derived by means of the Feynman- Vernon 
influence functional approach [20, 21, 24]. This approach 
allows exploring any regime, low or high temperature, 
strong or weak damping, etc. (For details on this approach 
see Refs. [20, 21, 24, 25] or Ref. [26].) At this point, we 
assume that initially the density matrix of S and TB 




FIG. 1. The system S, initially described by ps(0), is coupled 
to the thermal bath TB at temperature Ttb and frequency 
cutoff JItb- 7tb denotes the strength of the cooling. 



factorizes, i.e., p(0) — ps(0) ® Ptb(0) and additionally 
that each mode of the bath is at thermal equilibrium 
at the same temperature Ttb [20, 21, 24]. Within this 
framework, the time evolution of the density matrix is 
given by 

x J(q';,ql,t,q' + ,q'_,0)(q' + \p S (0)\qL), 

where the q label denotes the system coordinate rep- 
resentation, J(q+,q'_!_,t,q' + ,q'_,0) is the influence func- 
tional, which is given in terms of a path integral after 
tracing out the environmental degrees of freedom (see 
Refs. [20, 21, 24, 25] or Ref. [26]). This solution does not 
provide direct insight into the dynamics in the eigenba- 
sis of Hs- However, the solution given in Refs. [20, 21] 
can be analytically transformed into the eigenbasis of 
Hg (denoted by {\n}} with eigenvalues {E n }) giving the 
following expression: 

(n\p S (t)\m) = ^ Jnm;vv{t){v\fa{Q)\v) 

" (3) 

where J nm ;uu(t) is the influence functional in the energy 
basis representation, i.e., 

■W„(*) = Jdq'ldCdq' + dq'_J(q'l,C,t,q' + ,q'_,0) 

x{n\q'i){C\m){q' + \v)W_) 

The result provides a remarkable linear map between the 
initial system state and the final system state that (see 
below) hold in the case of initial system bath correlations. 

Of particular interest is the interplay between the diag- 
onal elements of the density matrix (state populations) 
and the off-diagonal elements (coherences). In absence 
of coupling to the environment the kernel J n m;v^{t) re- 
duces to J nm;vfl (t) = e-'^ Em - En)t 'H nv 8 mfl . The role of 
the 8 nv 5mp is twofold: it prevents the transfer of initial 
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population from (i/|/3g(0)|f) to (n\ps(t)\n) and addition- 
ally, it eliminates any effect of initial off-diagonal elements 
on the populations (i.e., the diagonal elements). By con- 
trast, in the presence of coupling to the environment, the 
overlap of the system eigenstates generates "new routes" 
for affecting populations, provided by terms of the type 
Jnn;w{t)- These terms transfer the initial population 
(v\ps{G)\v) to (n\ps(t)\n) at time t. Additionally, any 
initial off-diagonal elements contribute to the time de- 
pendence of the populations through the nonzero J n n;i/fj. 
terms. 

For the system-bath bilinear coupling in Fig. 1, the 
influence and nature of the bath is determined by a 
spectral density J(w) [20, 21, 24], which can be ex- 
pressed in terms of the bath modes parameters: J{uj) = 

c 2 

tt V* n — 1 — S(u! — u>j). In this case, we assume an non- 
Markovian Ohmic spectral density 



foj ^ TB = 8.2724 -- ^0/3^=1.0341 ■■- fiw /3 TB =0.5170 



J (to) = TO7 TB o;il| B / (fi| B 



(5) 



where 7tb is the strength coupling constant to the 
thermal bath and Htb is a frequency cutoff. This 
spectral density generates the following damping kernel 

"»7(*) = E, cos 



( Wj t) = 2/^^icosM) = 



7tb^tb exp(— OtbI^D) This kernel is responsible for the 
relaxation process and describes, roughly, the rate at 
which energy is transferred to the bath. In the limit 
when the cutoff frequency JItb tends to infinity, 7(4) — >• 
27tb<5(£), which corresponds to Markovian Ohmic dissi- 
pation. 

If initially ps(0) = |n)(n|, i.e., there is no initial 
coherence, the second sum in Eq. (3) vanishes and the 
time evolution of the diagonal terms depends exclusively 
on their initial values. For example, for the particular 
case when ps(0) = |0)(0|, the subsequent time evolution 
of the diagonal elements of the density matrix can be 
expressed simply by 



(n\p s (t)\n) = Jnn-oo(t)- 



(6) 



In Fig. 2, we have depicted J n n-oo(t) for some values of n. 
By contrast to the unitary case where J n n;oa(t) = S n oS n o, 
here we have the possibility of populating different en- 
ergy levels. As expected for a thermalizing system, the 
amount of population transfer increases with increasing 
temperature (see Fig. 2). So, it is clear that the role of 
Jnn-.oo(t) is transferring of population from (0|ps|0) to 
other eigenstates of S. This is a natural consequence of 
the fact that the spectrum of open quantum systems is 
broadened [20, 21, 24, 25], so excitation of single energy 
levels is not possible because of the overlap. These envi- 
ronmentally induced "new routes" for population transfer 
from different eigenstates have been extensively exploited 
in the context of biological systems [27, 28] . Here we can 
see that they appear mediating thermal activation, with 
the process being incoherent in nature. 

Remarkably, see Fig. 3, this population transfer is ac- 




FIG. 2. Time evolution of Joo;Oo(i) (blue curves), Jn ; oo(i) 
(green curves) and J22 ; oo(i) (red curves). Results are for 
7tb = O.Iojo, ^tb = 20ujo and hujo/ksTTB = 8.2724 (con- 
tinuous curves), fo^o/fcBTrB = 1.0341 (dashed curves) and 
hido/ksTTB ~ 0.5179 (dot-dashed curves). 



companicd by the generation of off-diagonal terms of 
the density matrix, (n\ps(t)\m) — J nm -oo(t). By con- 
trast to the diagonal terms, the asymptotic value of these 
off-diagonal terms are seen to increase with decreasing 
temperature and hence they become relevant at low tem- 
perature. For an harmonic bath, low temperature is 
accompanied by non-Markovian decoherence dynamics of 
the system. More importantly, the longer time decay seen 
in Fig. 3 is not related to decaying coherences because 
initially there was no coherence in the density matrix. 
Rather, they reflect the fact that in the effective basis, 
the populations are changing to reach the thermal state. 
It is worth noticing that these off-diagonal terms cannot 
be related to the existence of a coherent superposition 
of states, but to the overlapping of energy eigenstates 
induced to the the incoherent effect of the bath. 

Note, significantly, that terms like J nm -oo shown in 
Fig. 3 do not go to zero at long times. This implies that 
the associated [See Eq. (3)] (n\ps(t)\m) matrix element 
assumes a long-time constant nonzero value. 

In order to understand the nature of these off-diagonal 
terms, note that the density matrix is a double sided 
object, which means that we can interpret J nm ;Oo(t) as 
the influence functional associated with transitions from 
|0) — > 1 71) and (to I 4— (0| mediated by thermal activation 
(similar terms arise in optical nonlinear response [29]). In 
Fig. 4, we have schematically represented this situation: 
black lines depict the "symmetric" transitions, which are 
associated with transitions of the type |0) — > \n) and 
(n\ 4— (0|; colored dashed lines depict the "asymmetric" 
transitions, which are associated with transitions of the 
type |0) — > \n) and (to| <— (0|. Although these "asym- 
metric" transitions are less probable, they are enhanced 
at low temperature, fiu?o/3TB <C 1 and |S.7tb/3tb *C 1, as 
seen in Fig. 3. 

The fact that these off-diagonal terms survive at equilib- 
rium, points out the possibility of deviations from " canon- 
ical typicality" [30] (i.e., diagonal (n\ps(t)\n) at equilib- 
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fiw /3 TB = 1.0341 
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the effective frequency 




3.0 6.0 9.0 12.0 15.0 u t 

FIG. 3. Time evolution of Jo2 ; oo(i) (upper panel) and Ji3 ; oo(£) 
(lower panel). Results are for 7tb = O.Icjq, Otb = 20cjo and 
&Jo/fcB2TB = 8.2724 (continuous blue curves), fioJo/fceTrB = 
1.0341 (dashed green curves) and fejo/feTTB = 0.5179 (dot- 
dashed red curves). 



rium) at low temperature [31]. From our discussion, we 
can note that these off-diagonal terms are not attributable 
to interfering processes, but rather are part of thermal ac- 
tivation induced by the coupling to the bath, so they can 
take constant nonzero values at equilibrium rather than 
just vanishing. Based on this description, we should not 
call them coherences, but just static off-diagonal terms. 
However, these have been called stationary coherences, as 
noted in Sec. II and are simply static manifestations at 
equilibrium, of the system-bath coupling. 



.,•,.,1 = 2(/ l /?T B r 1 arccoth ( -yW) 



with (q 2 ) and (p 2 ) being the equilibrium variances [21, 32]. 
This allows expressing pp as [32] 



= Z 



/3tb 

n=0 



exp(-E n/3 PT B )\np)(np\, 



(7) 



with E„ 



Hui c ff(nj3 + h) being the eigenvalues and \np) 
the eigenstates of the effective Hamiltonian H c g. At high 
temperature, /iwo/^tb "C 1 and |S.7tb/3tb <C 1, m e ff and 
oj e g approaches their bare values m and ujq, respectively, 
and pp tends to the canonical distribution [32, 33] (with 
no off-diagonal terms in the system energy eigenbasis). 
At low temperatures, HluoPtb ^S> 1 and t^7tb/3tb ^> 1, 
they undergo to strong deviation due to damping [32, 33]. 

Since we are interested in natural processes, we next 
consider excitation of the equilibrated system S by (i) a 
second thermal bath TB' of the same nature as TB but at 
a different temperature and different coupling constant, 
and (ii) by blackbody radiation, denoted BB. For example, 
in the particular case of biological processes, e.g. in 
electronic energy transfer in photosynthetic complexes 
[1, 9], the coupling to the environment is strong [34] and 
the effective temperature is low[15, 16]. So, a question 
immediately follows: do these off-diagonal terms play any 
relevant role in any subsequent dynamics of these systems? 
More specifically, in what way do they contribute to any 
dynamics arising from subsequent perturbation? 
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FIG. 4. Symmetric and asymmetric transitions from the 
ground states |0)(0| to higher states during thermalization. 



IV. DYNAMICS IN THE PRESENCE OF A 
SECOND THERMAL BATH 

Equilibration of a system at temperature Ttb, with 
another bath at temperature Ttb' is of general interest. 
For example, for biological systems, thermal activation of 
biological processes by temperature changes can be found, 
e.g., in the context of transport of Ca ions through mem- 
branes [8]. To examine such processes, after thermalizing 
with TB, we couple the oscillator to a second dissipative 
environment TB'. The Hamiltonian of the system can 
now be written as 



H = H s + H- 



TB 



TB' 



■Hi 



ST 



Hi 



SB' 



(8) 



At equilibrium, the magnitude of the off-diagonal 
terms in the system density matrix can be evaluated 
from the equilibrium density matrix [32], i.e., pp = 

Zg^TrTB exp ^~/3 TB (Hs + H ST + # T b)) , where Zp TB 
is a normalization factor and Tttb denotes the trace 
over the bath. For this particular case we can introduce 
the effective Hamiltonian H c ff = ^ - p 2 + ^m c sLd 2 s q 2 , 



where the effective mass m e ff 



and 



where -tTtb' is the Hamiltonian describing the second 
thermal bath TB', while Hsb' describes the interaction of 
the system with TB'. As in the previous case, we choose 
Hbb' as being composed of a set of harmonic oscillators, 



H 



oo,2 



1_( 



q 3 ,k 



Cj ,kQx 



Note that the coupling is of the system S to the heat bath 
TB and of the system to the heat bath TB', as shown in 
Fig. 5. There is no direct coupling between TB and TB'. 



G 




FIG. 5. After thermalizing with TB, the system S is put in 
contact with a second thermal bath TB' (or BB, as in the 
plot) at different temperature and different coupling constant. 
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The resultant evolution of the system density-matrix 
can also be analytically obtained by using the influence 
functional approach [26]. The evolution is of the form in 
Eq. 3, but with ps(0) = pp and the J„ m;l , M (i)'s containing 
information about the effect of the initial correlation on 
the subsequent dynamics. 

In our approach, we have exact analytical access to 
every contribution to the dynamics. By contrast, for 
example, when applying the secular approximation as 
in Ref. [6], the contribution from the off-diagonal terms 
v 7^ p is ignored. 

It is important to note that in this case the initial 
condition is p(0) = ps+tb(0) ® /5 T b'(0)j where />s+tb(0) 
is the equilibrium density operator of (S+TB), and ptb' (0) 
denotes the thermal density operator of the second bath 
at temperature Ttb'- We describe the effect of the baths 
using, for both TB and TB', the functional form given 
in Eq. (5). In the absence of TB', the overall (S+TB) 
is time-independent, as it should be thermal equilibrium 
[21, 35, 36]. 

During this second thermalization step, the off-diagonal 
elements of the system density matrix that were generated 
during the first thermalization step affect the population 
dynamics (the diagonal terms) because they now enter in 
the initial density matrix [see Eq. (3)]. 

In the upper panel of Fig. 6, we present the time evo- 
lution of (0\ps(t)\2) and (l\ps(t)\3) for the conditions 
indicated in the figure caption. The off-diagonal terms 
are seen to reach the same order of magnitude as the 
populations (not shown), i.e. 10 _1 . Time evolution of, 
e.g., (0|/3g(t)|2) is affected by terms of the type J^.w and 
Jq2,uh- In particular, population transfer is assisted by 
Jq2,w and decay of the initial coherences is mediated by 
Jq2,vh- In the central and lower panels of Fig. 6 we have 
depicted Jq2,w and Jq2,vh, respectively, for various values 
of v and p. There we can see that while Jq2,w reach an 
asymptotic finite value, Jq2,v^ goes to zero as time evolves. 
Hence, it is clear that terms like J$2,vn are qualitatively 
different from Jq2,w type terms. The former are related 
to the decay of quantum coherences while the latter are 
related to population transfer during the thermalization. 

Figure 7 shows the time evolution of the ground state 




FIG. 6. Upper panel: exact time evolution of (0|j5s|2) and 
<l|ps|3> with riw /k B T T B = 8.2724, 7 tb = 0.1u) , fi TB < = 
20ojo, Ttb' = 2Ttb, 7tb = 27tb and f2 TB < = 2£!tb • Lower 
panels: time evolution of Jq2;v^ for some values of v and p. 



population (0|/Sg(£)|0). The exact evolution is depicted by 
the continuous blue curve, the evolution disregarding any 
correlation with TB (i.e., using the canonical distribution 
for ps(0)) is shown by the dot-dashed red curve and 
the evolution neglecting the off-diagonal terms (i.e. the 
secular approximation) in Eq. (3) is shown by the dashed 
green curve. The lower panels show some of the non- 
vanishing Joo;vfi(t) elements. The fact that the dynamics 
in these three cases is practically the same shows that, 
although off-diagonal terms are present in Fig. 6, they 
need not imply a significant contribution to the dynamics 
between eigenstates. That is, the presence of off-diagonal 
terms in the density matrix of open quantum systems in 
the system eigenstate representation is not necessarily an 
indicator of coherent effects. 

It is worth noting that off-diagonal terms generated by 
the presence of the second bath cannot manipulate the 
dynamics of the populations by themselves. This can be 
seen when projecting in the basis in which pp is diagonal, 
i.e. {n^ps^np) =Y JVft Jn f3 n f3 -,u f3 u f3 {t){v \p fj \v p ). This 

means that although off-diagonal terms are generated (see 
Fig. 6), they do not participate in the evolution of the 
populations. Only the combined action of the thermal 
baths, applied sequentially, can lead to the possibility of 
altering the dynamics of the populations. Connecting the 
system to the baths, either separately or simultaneously 
only leads to excitation of the system. 
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FIG. 7. Upper panel: exact (continuous blue curve), secular- 
approximated (dashed green curve) time evolution of (0|ps|0). 
The evolution using the canonical distribution as the initial 
state is depicted by the dot-dashed red curve. Lower panels: 
time evolution Joo;i^ for some values of v and fj,. Parameters 
as in 6. 
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where e/c is the coupling constant to the radiation, a^ s 
and a k are the annihilation and creation operators of 
the field mode of momentum k and polarization s. The 
vector potential is given by 



^* = E 

k,s 



he 


2 r 


W_ 





/feak, s e k!S q 



A«L e k, s • q > ( 10 ) 



where e is the polarization vector, V is the volume of 
the cavity and is the electron form-factor (Fourier 
transform of the charge distribution) which incorporates 
the electron structure [22]. 

By defining = Aire 2 /| / '(oj-^V) and ak, s = 

(^k^k9k.s + ipk,s)/ V^Wkfa^k, and applying the Power- 
Zienau's transformation (see Ref. [26] for details) we can 
rewrite Eq. (9) as 
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V. DYNAMICS INDUCED BY BLACKBODY 
RADIATION 



Consider now the case where the second thermal bath 
TB' is replaced by blackbody radiation (BB). Here, a 
charged harmonic oscillator is immersed in a dissipativc 
environment TB and coupled via dipole approximation to 
blackbody radiation BB. This provides a generic model 
for a wide variety of objects such as atoms, ions, electrons, 
molecules in equilibrium, subjected to blackbody irradi- 
ation. Qualitative results for this case are particularly 
relevant for natural light incident on biomolecules. 

The Hamiltonian of the total system is now of the 
form of Eq. (8) where Hb is the Hamiltonian describing 
the radiation field and Hsb' describes the interaction of 
the system with the blackbody radiation. In the dipole 
approximation, the Hamiltonian for an oscillator inter- 
acting with a radiation field and coupled linearly to its 



where the oscillator is seen to be coupled to the momentum 
coordinate pk,s- Since E x = —dA x /dt and based on the 
definition of the annihilation and creation operators (see 
above), one can show that this momentum coupling is 
equivalent to the oscillator being coupled to the electric 
field of the radiation with 

^ i 

Ex = ! E ( ~y ) (A* a k, s e kiS • q + /fc4, s e L • q) • 

k,s ^ ' 

In order to describe the action of the blackbody radia- 
tion, we need to consider, as in the previous case, that the 
modes in the cavity are thermally populated (for more 
details see Ref. [26]), i.e., each mode is characterized by an 
incoherent density operator. This feature introduces the 
incoherent nature of the radiation considered in this work. 
Based on this description, we expect that the blackbody 
radiation generates an incoherent evolution in contrast 
with the coherent evolution induced by laser pulses. 

Interestingly, the thermal fluctuations generated by 
the blackbody radiation are characterized by a two-point 
electric-field correlation function (E x (t")E x (t')) that is 
not (5-correlated. This reveals the intrinsic non-Markovian 



character of the radiation from both a statistical viewpoint 
[22, 26, 37] and from the optics point of view [38]. The 
correlation time t bb of the randomly fluctuating electric 
field can be calculated from the thermal time t bb — 
h/k^TsB [38] and is expressed roughly as t bb ~ t bb . 
For blackbody radiation at Tbb = 300 K, t bb ~ 160 fs 
whereas for sunlight, Tbb = 5900 K, t bb ~ 8 fs and 
for moonlight, T BB = 4100 K, r£ B - 12 fs. Thus, for 
processes taking place on the order of, e.g., 1 ps (such 
as electronic energy transfer in photosynthetic complexes 
[1, 9]) this coherence time is very short. Hence, under 
illumination by sunlight, the perturbation is effectively 
CW and incoherent. 

From an open-quantum-system perspective, the influ- 
ence of the blackbody radiation on the system is condensed 
in the spectral density [22, 37] 

j BB {uj) = MT BB uj 3 n 2 BB /(n 2 BB + Lj 2 ), (12) 

where M = m + A/tbb^bb is the renormalized 
mass, t bb = 2e 2 /3M c c 3 ~ 6.24 x 10~ 24 s and Q BB 
is a frequency cutoff [22]. This spectral density 
generates the following dissipative kernel 7bb(£) — 
tbb^Ib [2<J(t) — ^bb exp(— ObbI*!)] ■ Note that there is 
a fundamental limitation to the use of Eq. (12). That 
is, in the limit J7bb oo, we get the surprising result, 
7bb(£) = 0, i.e. no relaxation. This corresponds to 
the point-electron limit [//. = il BB /(f2 BB + lj 2 .) = 1 in 
Eq. (10)] and is unphysical because even for the electron, 
JIbb remains finite, although large. There is a natural 
upper limit given by [39, 40] J7bb = Tbb -1 , which cor- 
responds to 2/3 of the time for photon to traverse the 
classical electron radius (r°j = 2.818 x 10~ 15 m). Beyond 
this natural limit, causality is violated [39] and the bare 
mass m takes negatives values [39]. 

The exact analytic expression for the influence func- 
tional for this case is derived in Ref. [26]. There we show 
that the equations of motion are driven by a transient 
term that is dependent on the initial conditions and pro- 
portional to 7bb- This term is absent in the former case 
above, and is a consequence of the coupling to the mo- 
mentum of the modes rather than coupling through the 
coordinate of the modes [cf. Eq. (11)]. So, in addition to 
the turn-on effect, present in the former case, in the case 
of incoherent excitation by blackbody radiation, we also 
have a driven transient term. 

In order to gain insight into the strength of this transient 
term and of the blackbody radiation, we examine results 
in the limit when Qbb ~~ * t bb , and where the effect of the 
radiation can be estimated by a constant damping kernel 
given by 7bb = ^o t bb [22]. Assuming a typical value 
of electronic molecular structure (e.g. a carbon-carbon 
bound) uj — 3x 10 14 Hz, we have that 7bb = 1-8 X 10 -9 cj - 
In this case the excitation due to the radiation is too 
weak to compete with the incoherent effect of TB and, in 
absence of TB, the thermalization would take far too long 
(7 BB = 1.85 fxs) to be appreciable on, e.g., a picosecond 
time-scale. In order to see some sort of appreciable effect 



on the picoseconds time scale would require a cutoff on 
the order of J7bb = 5 x 10~ 6 t bb (7 bb = 9 ps). 

As an example of the dynamics, Fig. 8 shows the time 
evolution of the ground state and some propagating el- 
ements Joo,vfi(t) and Fig. 9 provides the corresponding 
results for the first excited state. In order to see some sort 
of appreciable effect, we have used here the artificial value 
of £!bb = 5 x 10~ 6 t bb and additionally strong coupling 
to the environment 7xb = 10 _1 wo- With these param- 
eters, we have that in Fig. 8 and Fig. 9, T T b = 277 K 
while Tbb = 5900 K. For this frequency, the time interval 
depicted in Fig. 8 and Fig. 9 is ps 128 fs. Here, the exact 
evolution of (0|ps (*) |0) and (l\ps(t)\l) are depicted by a 
continuous blue line, the time evolution neglecting the off- 
diagonal (stationary coherent) terms initially generated 
by the presence of the TB is depicted using the dashed 
green line, and the evolution using the canonical distri- 
bution as the initial state is shown in red. That is, as 
above, the results in red take the canonical distribution 
as the initial state whereas the dashed green calculation 
uses the thermal state (Eq. 6), but disregards the static 
coherences. 

In the lower panels we present some of the non- vanishing 
Joo-,vij,{t) and Ju- }V ^{t) elements. It is evident that the 
effect of the radiation is so weak that it generates nei- 
ther off-diagonal terms nor significant changes in the 
populations dynamics. That is, even with the choice of 
excessively aggressive parameters, the effect of the radi- 
ation is negligible. Rather, the time evolution can be 
seen as being caused by the turn-on of the interaction 
followed by the subsequent relaxation of the system to 
equilibrium, dominated by the interaction with TB. The 
weak oscillations in Fig. 8 and Fig. 9, can be attributed to 
the transient driven term, which is strongly determined by 
the frequency cutoff J7bb ■ For different sets of parameters 
these oscillations might well be absent. 

The effect of switching the interaction on t = and 
the time-dependent driving transient term can be also 
appreciated by the jump of the J„ m:t/A1 -terms. That is, 
in Fig. 8, Joo:00 is seen to jump from unity at t = to a 
different value, and Joo;ii and Joo;^ j um P from zero to 
a finite value. The same applies for Jn ; oo and J\i- VIJt in 
Fig. 9. 

From Fig. 8 we can see oscillations in the population 
of the ground state, these oscillations are also present in, 
e.g., the population of the first excited state [see Fig. 9]. 
Since, in our "secular approximation" these oscillations 
are removed, according to us, this would imply that the 
excitation is coherent, but very short lived. 

At this juncture, it is illustrative to comment on the 
difference between these results and those resulting when 
the equilibrated system (S+TB) is excited by a coherent 
source where it is assumed that it is the system S that 
interacts with the radiation. Given this coherent excita- 
tion (e.g. pulsed transform limited laser), the dynamical 
features change completely: the absorption of one photon 
from a coherent pulse creates a coherent superposition of 
energy eigenstates, and hence a time evolving state (cf. 
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Ref. [7]). It is worth mentioning that in absence of TB, if 
the initial state is a pure state, then this superposition 
will be described by a pure state. By contrast, under 
the same circumstance (absence of TB) , the blackbody 
radiation will create a incoherent superposition of states 
leading to a mixed state. In the presence of TB, the sta- 
tionary coherent terms created during the thermalization 
with TB will allow us to enhance the coherent control 
over the populations by transferring the coherence of the 
pulse into the populations [41]. 

As noted above, the parameters in the computation in 
Fig. 8 and Fig. 9 are artificial, and designed to show some 
effects due to BB. Based on naive classical arguments, this 
value of fieB would be equivalent to the inverse of the time 
needed by light to travel over around 1 A. In the context 
of electronic energy transfer in photosynthetic complexes, 
where natural frequencies luq = 6.63 x 10 12 Hz occur, 
then the coupling constant will be sa 7bb = 8 x 10 _6 u;o 
(7bb = 19 ns). In this case the dynamics induced would 
be far less noticeable than the one shown in Fig. 8 and 
Fig. 9. That is, since 7bb is smaller, then the effect on 
the relaxation will be weaker as well. Indeed, the final 
stationary population reached after equilibration will be 
far closer to that reached during the first thermalization 
step. Once again the response is to the turn-on of the 
field with rapid equilibration following mediated by the 
driving initial-condition-dependent-transient-term. 




body radiation is correctly represented as a thermal bath 
[see Eq. (12)]. Any time evolution observed here can 
be understood as the ultrafast transients as the system 
thermalizes, mediated by the driving initial-condition- 
dependent-transient-term. This term is present given the 
radiation field case, but is absent when excitation results 
from a second thermal bath coupled through the position 
of the bath modes. 
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FIG. 9. Parameters as in 8. 
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FIG. 8. Parameters for coupling to TB as in 7, with 
Huo/ksTsB = 0.3884, and Q B b = 8.3 x 10 s w . 

In summary, irradiating with blackbody radiation can- 
not generate a coherent dynamical result, since the black- 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

We have considered the natural process where a system, 
originally coupled to a thermal bath, is subsequently per- 
turbed by either another thermal bath or by blackbody 
radiation. We have shown that the first step of thermal- 
ization generates off-diagonal stationary coherences which, 
in principle, could affect the dynamics of the second per- 
turbative step. In particular, these off-diagonal stationary 
coherences could affect the populations of states during 
the second step. This process and is enhanced in the 
low temperature regime, and plays a fundamental role 
in explaining the origin of one photon phase control in 
molecular systems [41]. However, although our formal ap- 
proach clearly identifies the stationary off-diagonal system 
matrix elements as participants in the subsequent time 
evolution of the populations of a generic open quantum 
system, when these off-diagonal terms are generated in 
practice, i.e. in natural environments or by blackbody 
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irradiation, their contribution is negligible. This leads 
us to the conclusion that under natural conditions, the 
off-diagonal elements generated by thermal baths do not 
play any relevant role in the dynamics and that incoherent 
excitation of an open quantum system leads to dynamics 
free of coherent time evolution after an initially short 
transient time interval. 
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